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The quantum phase transition in an atom-molecule conversion system with atomic hopping be- 
tween different hyperfine states is studied. In mean field approximation, we give the phase diagram 
whose phase boundary only depends on the atomic hopping strength and the atom-molecule energy 
detuning but not on the atomic interaction. Such a phase boundary is further confirmed by the 
fidelity of the ground state and the energy gap between the first-excited state and the ground one. 
In comparison to mean field approximation, we also study the quantum phase transition in full 
quantum method, where the phase boundary can be affected by the particle number of the system. 
Whereas, with the help of finite-size scaling behaviors of energy gap, fidelity susceptibility and the 
first-order derivative of entanglement entropy, we show that one can obtain the same phase bound- 
ary by the MFA and full quantum methods in the limit oi N ^ oo. Additionally, our results show 
that the quantum phase transition can happen at the critical value of the atomic hopping strength 
even if the atom-molecule energy detuning is fixed on a certain value, which provides one a new way 
to control the quantum phase transition. 

PACS numbers: 03.75.Hli, 05.30.Rt, 05.30.Jp 



I. INTRODUCTION 

The quantum phase transition (QPT) describes an 
abrupt change in the ground state of a many-body sys- 
tem as some system parameters going across a critical 
point (at zero temperature). Quantum phase transitions 
(QPTs) in the systems of the quantum Hall, supercon- 
ductor and ultracold atoms have been studied extensively 
[ll-Q. Ultracold atomic systems, especially Bose-Einstein 
condensates (BECs), provide us a good platform to study 
the QPT. The experiment on ultracold atoms in an op- 
tical lattice by Bloch et al. has given a good example of 
the QPT from a superfluid to a Mott insulator (SF-MI) 
Then the SF-MI transition in the ultracold atomic 
systems was discussed widely [1, 0]. Besides the SF- 
MI transition, the transition from non-entangled to en- 
tangled states in two- mode BECs 0], from degenerate 
and non-degenerate ground states in the extended boson 
Josephson-j unction model Q, and from a pure molecule 
state to a mixed atom-molecule one in an atom-molecule 
model have been investigated. Note that in the 

above QPTs, one can control the QPT by changing dif- 
ferent system parameters such as atom-pair tunneling 
strength and energy detuning between the atomic and 
the molecular states. 

Molecular BECs are versatile in physical studies be- 
cause, in comparison to atomic BECs, they have more 
freedoms to be controlled. For example, the ultracold 
polar molecules have been used to study the ultracold 
chemistry, quantum many-body physics and quantum in- 
formation science. The polar molecules with dipole mo- 
ment may be either heteronuclear molecules or homonu- 
clear molecules. The heteronuclear molecules with large 
electric-dipole moment have been investigated both in 
theory and experiment widely [l^ - [l7| . The homonuclear 



molecules with large dipole moment, such as Rydberg 
molecules, have been produced by T. Pfau et al. in ex- 
periments P and studied widely jUill. One Ryd- 
berg molecule consists of a single kind of atoms in differ- 
ent states, and the atoms can jump between the ground 
state and Rydberg states, which provides probability to 
study the QPT in atom-molecule conversion system with 
atomic hopping. Although the quantum phase transition 
of the ground state has been studied in the atom-molecule 
system [9l-[l3| , the atoms in that system were in the same 
state and the hopping between different hyperfine atomic 
states was not considered. So it is need to consider the 
atom-molecule conversion system with atomic hopping 
and study the effect of the hopping strength between 
different hyperfine atomic states on the quantum phase 
transition. 

In this paper, we study the quantum phase transition 
of a atom-molecule conversion system with atomic hop- 
ping between different hyperfine atomic states in both 
mean field approximation (MFA) and full quantum meth- 
ods. It is interesting to find that the QPT can still ap- 
pear in the system as the hopping strength increasing 
even if the atom-molecule energy detuning is fixed on a 
certain value, which is different from Ref. where the 
QPT was induced by changing the energy detuning. In 
MFA, we show that the QPT exists in the thermody- 
namic limit (i.e., TV — > oo) and give the phase diagram 
whose phase boundary only depends on the atomic hop- 
ping strength and the atom-molecule energy detuning but 
not the atomic interaction. In full quantum method, we 
characterize the QPT with the help of energy gap, fidelity 
susceptibility and the first-order derivative of entangle- 
ment entropy, which give the same phase boundary. We 
also show that one can obtain the same phase boundary 
by the MFA and full quantum methods in the limit of 
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N ^ oo through studying the finite-size scahng behav- 
iors of energy gap, fidelity susceptibility and the first- 
order derivative of entanglement entropy, which further 
confirm the existence of the QPT in the atom-molecule 
conversion system with atomic hopping. 

In the next section, we give the model of the system 
and the general phase diagram in mean field approxi- 
mation. The energy gap and fidelity are also studied to 
characterize the QPT. In Sect. Hill we study the QPT 
from a pure molecule state to an atom-molecule mixed 
state in full quantum method. The energy gap, fidelity 
susceptibility, entanglement entropy, first-order deriva- 
tive of entanglement entropy, and their scaling behaviors 
are investigated. In the last section, we give a brief sum- 
mary. 



II. MODEL AND PHASE DIAGRAM 

We consider a three-component atom-to-molecule con- 
version system where the atoms can jump between two 
hyperfine atomic states. 

H = — J{a\a2 + dlai) + -^(alai — a^a2) 

+g'(b^aid2 + a\dlb) — ^(d\di + 0202 — b^b),{l) 

where cii and dl {i = 1,2) denote that annihilate and 
create a atom in the ith hyperfine atomic states, and 
b and 6^ denote that annihilate and create a molecule, 
respectively. Here J refers to the hopping strength be- 
tween the two atomic components, g' the atom-molecule 
coupling strength, Ua {U^) the strength of atomic intra- 
component (intercomponent) interaction, 6a the energy 
detuning between the two atomic components, and db the 
energy detuning between atomic and molecular compo- 
nents. 

In order to study the property of the system we con- 
sidered, we need to give the dynamical equations of the 
system in mean- field approximation (MFA). We know 
that with the help of Heisenberg motion equation for op- 
erators, one can easily give the evolution equations for 
operator di and b. In mean-field approximation, the av- 
erage value of the operators can be replaced by their av- 
erage values, i.e., a\ — {di)/^/N, a2 — {d2) /\fN and 
/3 = (5)/a//V- Then we can obtain the evolution equa- 
tions for ai , a2 and /3 with the help of the aforementioned 
evolution equations for di and b. Here ai and /3 satisfy 
the conservation law jaip -I- |a;2p + 2|/3p = 1 and N 
is the total particle number. For convenience to study 
the properties of fixed points, we assume ai — y/pie^^^ , 



"2 = ^/pi& and /3 = y/pb^^ > where pi + p2 + 2pb = I. 
The mean-field dynamical equations of the system can 



be written as 
d_ _ ^ J cos{2(j)a) - g^/p^ cos (j) t{ _ t^i ^ 

di"^"" ~ ^ ^{i-2pbf-z^ + 2" " + y ' 



dt ' 



2Jcos(20a) — 2gy/pbCOS( 

V/(1 - 2pb)2 - 



- U - U']il - 2pb) 



dt 
d 

di 



Pb 



g (l-2pb)2-z2 36b 
+ 2V 7b ^°^'^+-' 

: -2JV(l-2p6)2-z2sin(20,), 
: g^/ ((1 - 2/9b)2 - z^)pb sin0. 



(2) 



where 0a = (02 - di)/2, 4> = Bi + 62 - Ob, z = pi - 
P2. Here the renormalized parameters are U = NUa, 
U' = NUa and g = VNg'. The equations (0) satisfy 
the Hamiltonian canonical equations, i.e., ■^^'Pa — ^^f-^ 
= -ff' Tt<f>-^ and j,pb = where the 

Hamiltonian in MFA reads 



HmF = -Jv/(1 - 2/96)2 - z2 cos(20„) -I- -j-z 

+ ^iil-2pbr + z') + ^{{l-2pbr-z') 



+g^[{l - 2pbY - z-^]pb co&4>+ ^(3p6 - 1). 

(3) 

Since the fixed points correspond to the eigenvalues of 
the nonlinear system [l^l, we can study the static prop- 
erties of the system by solving the fixed points. The 



fixed points requires that = 0, ^0 = 0, J^z = 0, 

■^Pb — 0, which can give many fixed point solutions. In 
the following discussion, we only focus on the kind of 
fixed points with z = since they contain the ground 
and the first excited state. Such a kind of fixed points 
read as, 



Pb 



A± 4-652 



65 cos (j) 
z = = 0(7r) 



20a = 0(7r), 



(4) 



where A = 3(5h/2 -|- 2Jcos(2(/)a), is taken for = 0, 
and ' — 'is taken for = tt. The Equation @ is obtained 
in the assumption of 5i2 — and U ~ U + U' — 0. 
Note that for the case of f/ ^ 0, we can not give the 
analytical solution of the fixed point. Whereas, we find 
that the existence of U does not change the boundary 
between the pure molecule (PM) phase and the mixed 
atom-molecule (AM) phase. 

In Fig. [TJa), we plot the dependence of density of 
molecules on hopping strength for the ground and first 
excited states. In such a figure, the three quantities in the 
legend denote the values oiU, 20a and 0. The lines with 
{U, 0, tt) correspond to the ground states for different U, 
and the line with (0, tt, tt) corresponds to the first excited 
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FIG. 1. (Color online) (a) The density of molecules pb versus 
J of the fixed points with z — for [7 = marked with 
four cyan lines. The parameters are g = 1 and Si, — —2. U, 
2(j)a and <j) are indicated in the legend. The ground states 
for nonzero atomic interaction U = 2, 4, 6, 8 are also shown, 
(b) The phase diagram of the ground state in MFA in the 
parameter space of Sf, and J with unit of g. The pure molecule 
(PM) phase and the mixed atom-molecule (AM) phase are 
separated by a critical line (marked with blue solid line). We 
focus our thoughts on the red dotted line and the critical 
point {{5b = — 2(7, Jc = [3 — \/2\g/2) marked with a black 
filled circle). 



Fig. [2l we find that the energy gap between the ground 
and first excited state is zero for J < and nonzero for 
J > Jc for any interaction strength U. That implies that 
the transition from PM phase to AM phase happens on 
the point where the energy degeneracy between ground 
and first excited state is lifted [1]. We know that the 
fidelity for eigenstates can also be used to characterize 
the phase transition. Here the fidelity for the ground 
state is defined as Fmf = (4'gs(>/)|^'gs(>/ + ^J)), where 
|*Gs(>/)) and |^'gs(-'^ + SJ)) is the ground state of the 
system for different hopping strength J and J + SJ with 
dJ being a very small quantity. From Fig. [21 we also find 
that the fidelity has a sudden decrease from one at the 
critical point J = Jc, where the phase transition from 
a pure molecule phase to a mixed atom-molecule phase 
takes place. In a word, both the energy gap and the 
fidelity have fantastic phenomena at the critical point. 
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FIG. 2. (Color online) Atomic hopping strength (J) profiles of 
the energy gap between the first excited state and the ground 
state (A^mf) and of the fidelity of the ground state in MFA. 
The parameters are &b = —2(7 and JJ = O.OOlg. 



state for [/ = 0. Note that the_first excited state can not 
be influenced by the value of C/, so the line with (0, tt, tt) 
can also denote the first excited state for different inter- 
action strength U . From Fig. [TJa), we can find that the 
ground state is a pure molecule state when J < Jc, and 
it is a mixed atom- molecule state when J > 3^. Here the 
critical value Jc = {^ — \fi)gj2 which is dependent on the 
energy detuning (here (5{, = —2(7) but not on the inter- 
action strength U . Meanwhile, we give the general phase 
diagram of the ground state in the parameter space of (5^ 
and J in Fig. [T] (b). The boundary between PM phase 
and AM phase is 3(5fc/2 -I- 2 J = — v^.g with g = 1 (corre- 
sponding to the blue solid line in Fig. [1] (b)). Note that 
for the case of J = 0, the critical point (5j, = —\/sj9g) 
which agrees well with the result in ref. |l2| . 

In order to further confirm the phase boundary, we plot 
profiles of the energy gap between the first excited state 
and the ground state (Ai^Mp) and of the fidelity of the 
ground state for a fixed energy detuning db in MFA. From 



III. QUANTUM PHASE TRANSITION 

To get insight into the QPT, we study the system in full 
quantum method. For a finite particle number N, where 
one molecule is counted as two particles, the Hamilto- 
nian can be exactly diagonalized on the basis of the Fock 
state \n, N — n — 2m, m). Here n is the number of atoms 
for the first atomic component and m is the number of 
molecules. The dimension of the Fock basis is (A^/2-|-l)^. 
For convenience, the form of the Fock basis is signed as 
|n, m). Then the eigenstate of the system can be written 
as 

|*(J))= CnMJ)\n,m), (5) 

n+2m<N 

where Cn,m{J) are complex coefficients with parameter 
J. The eigenvalues and the eigenstates of the system as a 
function of J can be easily obtained based on the method 
of exact diagonalization of Hamiltonian 
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FIG. 3. (Color online) The total population of atoms {pa — 
1 — 2pi,) in the ground state as a function of the hopping 
strength J for different particle numbers. The red solid line 
is the MFA result {i.e. N oo). 



In Fig.[3l we plot the total population of atoms for the 
ground state of the system with finite number of particles 
and in MFA . From this figure, we can see that the QPT 
from a pure molecule state (the total population of the 
two atomic states is zero, i.e., pa = 1 — 2pb = 0) to 
a mixed atom-molecule state (the atomic population is 
nonzero, i.e., pa = 1 — 2yOf, > 0,) happens in the system 
as the atomic hopping strength reaches the critical value 
of Jat. Note that the critical value of Jjv is very close 
to Jc obtained in MFA, and Jn ^ Jc in the limit of 
iV — > oo. In the following discussion, we will further 
characterize the quantum phase transition with the help 
of the concepts of energy gap, fidelity susceptibility and 
entanglement entropy. Additionally, we will study the 
scaling behaviors of the system near the critical point Jc . 



A. Energy gap 

Now we are in the position to study the QPT with the 
help of energy gap which has been used to characterize 
the QPT in ref. [lll - ITsl . ITU . Here the energy gap is de- 
fined as the energy difference between the first excited 
state El and the ground state Eq, i.e., AE = Ei — Eq. 
In Fig.|l](a), the energy gap AE is shown as a function 
of J for different particle numbers. For a fixed energy de- 
tuning 6b, the avoided level-crossing between the ground 
and the first excited state appears near the critical hop- 
ping strength Jc which is the phase transition point given 
in mean-field approximation. From Fig. 2] (a), we can 
see that for a finite particle number N, the energy gap 
reaches its minimum value AEmin at the critical points 
Jn where the QPT from a pure molecule state to a mixed 
atom-molecule state happens. We find that the critical 
hopping strength becomes closer with the increase of the 
particle number N . 

To further characterize the finite-size effect present in 
Fig. m (a), we show the scaling behaviors of AEmin and 



AJ = Jm — Jc on the particle number N for different 
interaction strength in Fig. H] (b) and (c). In such two 
panels, the discrete points denotes the numerical results 
for finite particle numbers N , while the solid lines are 
the fitting functions which well show how the quantum 
results tend to MFA ones with the increase of particle 
numbers. From Fig. U (b) and (c), we find that both 
AJ and {AE)min converge to zero with different slopes 
for different U when N~^'^ and approach zero. So 
that both A J and [AE)min approach zero in the limit of 
N GO, i.e., AEmin converges to zero and Jn converges 
to Jc, which is agreed well with the results in MFA. Ad- 
ditionally, the critical exponents 5g and p, for different U 
with slight difference are shown in Table HI 



TABLE I. Finite-size scaling behaviors of the energy gap. 
Scaling exponents 5a of AJ ~ N~^'^ and p of AEmin ~ A'^"^ 
at various U obtained by sampling system size in A G 
[60,200]. 
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B. Fidelity susceptibility 

We know that the quantum fidelity can be used to 
characterize the quantum phase transition (ill . Iisl |23| - 
l25j , where the quantum fidelity is defined as the absolute 
value of the overlap between two ground states with an 
infinitesimal variation of the control parameter. In order 
to study the effect of hopping strength J on the QPT of 
ground state, the fidelity can be written as 

FiJ,dJ) = \{^o{J)\^oiJ + SJ))\, (6) 

where |'I'o(J)) and \'i>o{J + SJ)) are two ground states of 
the system with small parameter difference SJ. We find 
that the value of fidelity is dependent on the value of 6 J 
although there is a sudden drop of the fidelity value near 
the critical point Jc which is the phase transition point 
given in MFA. In order to make up for this deficiency 
(i.e., the value of fidelity is dependent on the value of SJ) , 
we make use of the concept of fidelity susceptibility [8l.l23l. 
|26| . In the first-order perturbation theory, we consider 
the hopping term as the perturbation term {i.e., H = 
Ho + JHj). Then the fidelity susceptibility is defined as 
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FIG. 4. (Color online) (a) Energy gap between the first excited state and the ground state versus the hopping strength J for 
different particle numbers, (b) The finite-size scaling behaviors of AJ — Jjv — Jc for different atomic interaction U. Jn is the 
parameter corresponding to AEmin- (c) The finite-size scaling behaviors of AEmin for different atomic interaction U. The 
parameter is Si, = —2g. 




J (unit of g) N-Sp (NT' x10"' 

FIG. 5. (Color online) (a) The fidelity susceptibility for the ground state as a function of the hopping strength J for different 
particle numbers, (b) The finite-size scaling behaviors of AJ = Jn — Jc for different U. Jn is the parameter corresponding to 
Xmax- (c) The finite-size scaling behaviors of Xmax for different U. The parameter is St — —2g. 



which does not depend on the value of 5J . It is also 
proved to be related to the correlation function which 
is used to show phase transition. 



The numerical results of the fidelity susceptibility ver- 
sus J for = 10, 50, 110, 190 are shown in Fig. E] (a). 
We can see that the fidelity susceptibility is about zero 
when J is far away from the critical point Jat, while it 
increases suddenly and reaches its maximum value at the 
critical point Jn which is dependent on the particle num- 
ber N . Meanwhile, we can find that the maximum values 
of the fidelity susceptibility become larger and the crit- 
ical point Jis! become closer to Jc as N increasing. In 
order to study the effect of particle number on the criti- 
cal behavior, we show the finite-size scaling behaviors of 
A J and {xmax)^^ for different U by power-law in Fig.[B] 
(b) and (c). 



We can find clearly that both A J and {xmax) ^ con- 
verge to zero for various U with different slopes when 
N~^'' and {N'^)~^ approach zero, respectively. In other 
words, A J is close to zero and Xmax is approximated to 
infinity when N oo. We show the critical exponents Sp 
and for different U without good convergency [1, [2^ , 
in Table ini 



TABLE II. Finite-size scaling behaviors of the fidelity sus- 
ceptibility. Scaling exponents Sf of AJ ~ N~^^ and u of 
Xmax ~ N'"^ at various U obtained by sampling system size 
in iV G [60,200]. 
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0.7005 


0.6781 


0.6716 


0.6682 


0.6659 
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1.336 


1.308 


1.297 


1.289 


1.281 



C. Entanglement entropy 

If the system can be viewed as a bipartite system, the 
entanglement entropy of its ground state is physically 
meaningful and has been used to characterize the quan- 
tum phase transition 0, [III [H [H, H, H^]. The von 
Neumann entropy, one typical entanglement entropy, of 
a bipartite system AB for a pure state \^) is defined as 

S = -TrA{pA log2 PA) = -Itb^Pb log2 Pb)-, (8) 

where Pa{b) = TrB(^) (|*)(^'|) is the reduced density ma- 
trix of the system with two subsystems A and B. In this 
paper, we consider the two atomic modes as subsystem 
A and the molecular mode as subsystem B. From the 
Schmidt decomposition of pure state, we know that the 
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FIG. 6. (Color online) The entanglement entropy (a) and its first-order derivative (b) of the ground state as a function of the 
hopping strength J for different particle numbers, (c) The finite-size scaling behaviors of AJ = Jjv — Jc for different U . Jjv is 
the parameter corresponding to dS / dJmax- (d) The finite-size scaling behaviors of dS/dJmax for different U. The parameter is 
5b = -2g. 



entanglement entropies calculated from the reduced den- 
sity operator of the subsystem A and B agree well with 
each other. 

With the help of the definition of entanglement en- 
tropy, we plot the numerical results of the entangle- 
ment entropy between the two subsystems for the exact 
ground state in Fig. [6] (a) with different particle num- 
bers N = 10, 50, 110, 190. In Fig. M (a), the maximum 
value of the entanglement entropy become larger as TV in- 
creasing. Whereas, unlike the behavior of the transverse 
Ising model (sil [33 ] , the entanglement entropy does not 
reach its maximum value at the critical point Jc even if 
N ^ oo. However, it is exciting to find that the sudden 
rise of the entanglement entropy takes place near the crit- 
ical point Jc = (3 — -\/2)/2. In order to describe the quan- 
tum phase transition, the first-order derivative of the en- 
tanglement entropy with respect to J is introduced. The 
numerical results of the first-order derivatives are shown 
in Fig. [S] (b). The maximum value {dS / dJ)max of first- 
order derivative of the entanglement entropy appears at 
the critical points Jm which become closer to Jc as N 
increasing. So as to find the dependence of {dS/dJ)max 
and A J on the particle number N, we show their scaling 
behavior for different U in Fig. [S] (c) and (d). We find 
that both AJ and {dS/dJmax)~^ converge to zero for 
various U with different slopes when N~^^ and {N'^)~-^ 
approach zero. So that A J is close to zero and dS/dJmax 
is in limit of infinity when N ^ oo. Meanwhile , the criti- 
cal exponents ^5 and lo for different U are shown in Table 



imi 

TABLE III. Finite-size scaling behaviors of the first-order 
derivative of the entanglement entropy. Scaling exponents 
Ss of AJ ~ N"^^ and w of {dS/dJ),^ax ~ N'^ at various U 
obtained by sampling system size in A*' G [60, 200]. 
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IV. SUMMARY 

In this paper, we have studied the phase transition in 
an atom-molecule conversion system where the atoms can 
jump between two atomic hyperfine states. In mean-field 
approximation, we have given the the phase diagram for 
the ground state, and shown that the phase boundary 
between pure molecule phase and mixed atom-molecule 
one is only dependent on the hopping strength J and 
energy detuning 5b but not dependent on the atomic in- 
teraction U. With the help of fixed points, we have stud- 
ied the fidelity for the ground state and the energy gap 
between the first-excited state and the ground one. We 
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have found that the ground state changes from degener- 
acy to non-degeneracy and the fidelity decreases suddenly 

at the phase boundary Jq = ~jSb ~ with g = 1, 
which implies that the energy gap and the fidelity can 
well characterize that phase transition. 

In comparison to mean-field approximation, we have 
investigate the quantum phase transition of the system 
in full quantum method. Taking the total population 
of atoms as the order parameter, we have shown that 
the QPT from a pure molecule phase to a mixed atom- 
molecule phase happens at the critical point Jjv for dif- 
ferent N with a fixed energy detuning Sh- Note that 
the value of Jn depends on both the value of Sb and 
the particle number N of the system. We have further 
characterized the QPT with the help of the energy gap 
between the first-excited state and the ground state, the 
fidelity susceptibility of the exact ground state, and the 
entanglement entropy and its first-order derivative be- 
tween the atomic subsystem and molecular subsystem, 
respectively, which give the same critical point Jn- We 



have also shown that the critical point Jjv approaches to 
Jc with the increase of particle number N through study- 
ing the finite-size scaling behaviors of the energy gap, 
fidelity susceptibility and first-order derivative of entan- 
glement entropy. So that in the limit of AT — >■ cxd, the 
MFA and full quantum methods can give the same phase 
boundary. Our results enrich the phenomena of QPT 
in multiple component system especially atom-molecule 
conversion system, and indicate that one can control the 
QPT of the atom-molecule conversion system by chang- 
ing the hopping strength between the two atomic hyper- 
fine states as well as the atom-molecule energy detuning. 
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